A central limit theorem for a localized version 

of the SK model 

Sergio de Carvalho Bezerra * Samy Tindel 

Institut Elie Cartan, Universite de Nancy 1 
BP 239, 54506-Vandoeuvre-les-Nancy, France 
[bezerra, tindel] Siecn . u-nancy . f r 

February 2, 2008 



Abstract 

In this note, we consider a SK (Sherrington-Kirkpatrick)-type model on Z*^ for 
d > 1, weighted by a function allowing to any single spin to interact with a small 
proportion of the other ones. In the thermodynamical limit, we investigate the 
equivalence of this model with the usual SK spin system, through the study of the 
fluctuations of the free energy. 
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1 Introduction 

This paper is concerned with a localized version of the Sherrington-Kirkpatrick model 
with external field, which can be described in the following way: for N,d > 1, our 
space of configurations will be S = Etv = {—1,1}'"^, where Cat is the finite lattice box 
Ctv = [—N] in Z*^. For a given configuration a G S^v, we will consider the Hamiltonian 

^""^ = ^ ^ ^ [~n) ^(*'^)''^'^^^ + ^ I] (1-1) 

where (3 stands for the inverse of the temperature of the system, = 2N + 1, 

is the notation for a pair of sites i,j E Cn (taken only once), {g{ij) '■ E Cn} is a 

family of IID standard centered Gaussian random variables, and h represents a constant 

*This author's research partially supported by CAPES. 
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positive external field, under which the spins tend to take the value +1. Our localization 
is represented by the function g, which can be thought of as a smooth frame, and which is 
only assumed to be defined on [—1, 1]'^ such that is of positive type, so that is non- 
negative and invariant by symmetry about the origin. We also assume q^ is a continuous 
function, including at its periodic boundary —1 = 1. The aim of our article is then to 
study the limit, when oo, of the Gibbs measure G'Af(o") defined on S at by 

„-£f]v(o-) 

Gjv(^) = , where Zn = Z^iP) = V e-^-^"^), 

and more specifically, we will concentrate on the so-called free energy of the system, 
defined by: 

pi/3) = lim E[pn{(3)] = a.s. - lim where = ^ log(Z^(/?)). (1.2) 



The model described by (jl.lj) can be considered as a finite range approximation of the 
mean field SK model, associated with the Hamiltonian 

-Hn (a) = 9{z,j)(^i'^j + hJ2^h 

for which a large amount of information is now available [H[ E] . It seems then natural to 
try to approximate the realistic spin glass system, on which we have very little rigorous 
knowledge (see however [7]), by our localized model ()1.1|) . capturing some of the geometry 
of the physical spin configuration, but still of a mean-field type in the limit N ^ oo. One 
could then hope to perform an expansion in in order to quantify the difference between 
the original SK model and our model (jl.lj) . 

In fact, this kind of idea is not new, and goes back at least, in the spin glass context, to 
A version of our model with /i = has been studied then in PU]; and more recently, 
the Kac limit of finite range spin glasses has been considered in j3 |2] . In these latter 
references, a slightly different point of view is adopted: the finite range model depends 
on a given parameter 7 > 0, (which would be 1/A^ in our setting), and this localization 
parameter is sent to after the thermodynamical limit in A^ is taken. It can be shown 
then, by some nice and soft interpolation arguments, that in the limit 7 — > 0, the free 
energy of the localized system is the same as the free energy of the SK model, for any 
value of the parameter /3 > 0. Notice that the results contained in 012] cannot be applied 
directly to our model, since in our case the limits 7 — > and N —>■ 00 are taken at the 
same time. However, some slight modifications of the computations contained in these 
papers would also show that our quantity Pn{/3) defined at (jl.2j) behaves like the free 
energy of the SK model for large A^. 

The goal of our paper is then, in a sense, more modest than E] , since we will only 
deal with the high temperature region of the model, i.e. small values of (3. On the other 
hand, our scope is to show that the equivalence between the SK model and our localized 
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model still holds, in the thermodynamical limit, for a second order expansion of the free 
energy, that is in the central limit theorem regime. More specifically, we will show the 
following limit result: let 70 be the L -norm of q in [— 1, 1]'^. For /3, /i > 0, let also s be 
the unique solution to the equation 

s = E [tanh2(/5^v^ + h)] , (1.3) 

and SK(/9, h) be the function 

SK(/?, h) = (3^1 - sf/A + log2 + E [log [cosh {/3zyG + h)]] , (1.4) 

which represents the free energy of the SK model in the high temperature region. Set also 

1 /2 

p{(3, h) = SK(7q j3, h). Then, under suitable conditions on g, we have 

(£)- lim iV^/2[p^(/5)-p(/3,/,)] = r, 

where y is a centered Gaussian random variable with variance r = r(/9, /i, 70). The 
announced equivalence, at the CLT level, between the SK model and the localized one, 
springs then from the fact that r(/?, /i, 1) is also the variance of the Gaussian random 
variable which shows up in the central limit theorem of the SK case (see IH E])- 

Let us say a few words about the method we have used in order to get our result: 
since we are in the high temperature regime, we are allowed to use a cavity type method 
in order to compute the limit of the overlap of the localized spin system. This yields then 
the limit of the free energy in a straightforward manner. It has also been shown in jH] 
that the stochastic calculus tools developed in P could be adapted to the case of spin 
glasses with external field. This induces a powerful method for obtaining central limit 
theorems for the free energy, and interestingly enough, in this context, a dynamical point 
of view gives some insight on a static stochastic problem. We will elaborate here on this 
method in order to treat the localized case, by taking advantage systematically of the 
Fourier decomposition of q. 

Our paper is divided as follows: at Sectional we compute the simple limit of the overlap 
function and of pn{P), recovering the results obtained in j2] for the high temperature 
regime. At Sectional we derive the announced central limit theorem thanks to stochastic 
calculus tools. 



2 Simple limit of the free energy 

Recall that we are dealing with the system induced by the Hamiltonian p.lj) . and let us 
define some additional notations about Gibbs averages: let / : — ^ M be a function of 
n configurations, with n > 1. Then we set 

^(/) = ^ E fi^\---,cTnexp[-J2HN{cT')], and u{f) = E [p{f)] . (2.1) 



1=1 
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In the sequel of the paper, we will also write du^ instead of dip/du for the derivative of a 
function tp with respect to a parameter tt, and Pn{P) = ^[pn{P)]- With these notations 
in hand, our strategy in order to get the limit of Pn{P) will follow the classical steps of 
the cavity procedure, namely: 

1. Find an expression for dppN{(3) in terms of an overlap-type function i?^'^. 



2. For a general function / : M, find a useful expression for 9„p„(/) along a 

suitable path defined for v G [0, 1], involving a Hamiltonian HN,v{cr)- 



3. Compute p^iR^'"^) inductively and deduce an expression for dppNiP)^ and then for 

We will start by the first of these steps, for which we will introduce a little more notation: 
first of all, we assume for the moment the following basic hypothesis on q: 

Hypothesis 2.1. The function q is continuous on [—1,1]'^. 

In this case, the function q^ : [—1, 1]"' R^ can be decomposed, as a function of 
L^([— 1, 1]"'), into a Fourier series of the form: 



q^{x) 



nrk-x 



with 7^ < oo and F = 7^ < 



00. 



(2.2) 



Set also R^'^ = iV""' Eiec^v ^i^i^ ^ ^ ^ ' 



R 



1,2 



X V ^ i7vi-k 19 /" ■» / 1 I-* 

>T 2^ e iv ajaf, lik^O, and i?, 



1,2 _ d1,2 




R^'^ - r, 



(2.3) 



iGCjv 



where r is a positive constant, whose exact value will be determined later on. Eventually, 
we will denote by the function q[-/N) defined on [— A^; A^]*^, and = Z'^\{0}. Then 
the following relation holds true: 



Proposition 2.2. For all P > 0, we have 



f 7.E [p{R]:')] - f 7oE [p{R^'')] + ^ 



Proof. Set B{(j) = e-^^^''\ Then the previous definitions and an elementary Gaussian 
integration by parts yield (see jHj): 



{i,j)eCN 



^ ^ qU^-j)cr }a}'yfa]Bia')B{a' 



72 



Y 4(^-j)E[p(a,V)afaJ)] 

(j,i)GCjv 
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Hence, using the decomposition (j2.2j) of g^, we obtain 



df3PN{P) 



JL 
J. 



■E 



E 



(ij)eCjv (i,i)eCjv Xkez-i 

2^ 



from which the desired expression is deduced easily. 



□ 



We are now ready to start the second step of the strategy mentioned above. 
2.1 The cavity method 

The cavity method will consist here in suppressing in a continuous way the interactions 
between a certain site m E Cn and the remaining spins. This will be done in the following 
way: set C"^ = {z G Cat; i ^ m}. We decompose, for all a G S^v, 

-Hn {(t) = -H"^ {pm) +am{h + Qm (Pm)) , 

where pm denotes the ordered spin values except for the m-th spin, that is, the (A^*^ ^ 1)- 
tuple (cTj; z G C]^) and we denote 



9m (Pn 

and we also use the notation 

-Hn (Pm) -- 



/3 



J\fd/2 



^ qN{m - i)g{m,i)(^i, 



This new Hamiltonian is similar but not identical to the Hamiltonian —H^-i on Tj^-i- 

The path we will build is now of the following form: consider a collection ); i G 

C*]^, V G [0, 1]} of independent standard Brownian motions, and {X{v)] v G [0, 1]} an in- 
dependent reversed time Brownian motion, all defined on the probability space JF, P). 
Notice that gii^rn) can be seen as the final value of the Brownian motion -Bj^m, and 

that X is the solution to a stochastic differential equation of the form 



X{v)=7]- [ ^^ds + W{v), ve[0,l] 
Jo i- — s 



(2.4) 
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where is a standard Gaussian random variable and W another Brownian motion, inde- 
pendent of the remainder of the randomness. Notice that, for notational sake, we have 

1 /2 

written Bi^^n instead of B^i^rn)- Set then r = 7o r and for v G [0, 1], define 
-HN,v{(r) = -Hn{p^) + ^^ ^ qN{m-t)a^amBi,m{v)+(3f^^^ ^ a^X{t) + h 

For / : — > M, denote also by Pvif) the associated Gibbs average, defined in a similar 
way to and i^m,v{f) = E[p^(/)]. 

Recall now the following elementary lemma from jH] : 

Lemma 2.3. For k > 1, let {Bf, I < k} be a collection of independent standard Brownian 
motions. Let also X be the solution to Ii2.4\ ), o-nd ip : R'^"'"^ M &e a function having 
at most exponential growth together with its first two derivatives. Then, for any v G [0, 1], 

E[^{B,{v),...,Bk{v),X{v))] 

k 

= E [^{v)] + ^ r E [dl^^^iB.is), . . . , B,is),Xi.s))] ds 

^-^0 1=1 



[dl^^^,^^^^{Br{s), . . . , B,{s),X{s)) 



ds. 



These preliminary tools yield the following differentiation rule: 



Proposition 2.4. For any f : M and v G [0, 1], the derivative of i'm,v{f) is given 

by 



keZ-i \l<l<l'<n 

E( ft?l^n+l I n+l\ I '^('^ + 1) ( ff?n+l,n+2n+ln+2\ \ ( o ^\ 

l^m,v l^fRk (^m(^m ) + ^ ^m,v [j Rk j I • (2-5) 



n 

n 

1=1 



Proof. This result stems from an easy application of Lemma 12. HI to the function 

G C]^; X,n{v)) = pmAf)- 

The computations of the second derivatives of ip are a matter of easy (though cumbersome) 
calculations, and are left to the reader for sake of conciseness. 

□ 



Now that the variations of Vm,v{,f) have been computed, we can proceed to get some 

k 



1 2 

bounds on the overlaps R ' 
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2.2 Bounds on the overlap 



Let us start with three lemmas whose proofs follow essentially that of the corresponding 
result in the cavity method for the standard Sherrington-Kirkpatrick model, and which 
will be combined with the explicit expression for i'm,v in Proposition 12.41 and a separate 
calculation of h'm,v for f = to obtain information on the actual expected overlaps, under 
= ^m,ii i-e. for V = 1. 

Lemma 2.5. There exist two positive constants Cn.q^p and c'^^p that depend only on n, q 
and P, and are uniformly bounded in [3 for P G [0, 1], such that if f is a positive function 
on then for all m E Cn , 



(/)< if) 



and 



k„.,.(/)-^„.,o(/)i<<,,^/3v/2(i/n 

Proof. See Propositions 2.4.6 and 2.4.7]. 



(2.6) 
(2.7) 

□ 



For V = 0, the expression of i^vif) can be simplified for a large class of functions / on 
Sat. 

Lemma 2.6. For fixed m G Stv, let f be a function on that does not depend on the 
values (7^, cr^, . . . , o"^. Then for any subset I of {1, ... ,n} we have 



/JI^L) =E[tanh(r)l^l]z/,„,o(/) 
i&i ) 



where Y is the Gaussian random variable defined as: 

Y = (3zy/^ + h, 

with a standard normal variable z. 
Proof. See jOl Lemma 2.4.4]. 



□ 



On the other hand, some symmetry properties for v = 1 yield the following kind of 
estimate: 



Lemma 2.7. Set 



5,= 5^7,z/„(|/?^'2|2), for t;e[0,l]. 



(2i 



Then 



S, = 70. {{R''^ - r) {alal -r)) + l-^Y.Y. ^^''^""""^ i^'^^^^) 
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Now, in order to exploit Lemma (2.61 we must modify the above expression for 6i by 
completing the following two tasks: 

(i) estimate the error made by replacing the arguments of Ui^Q by functions that are of 
the same form as those in Lemma I 



(ii) estimate the error made by replacing z/ by u^fi (or z/j o as appropriate). 

2.2.1 Task (i). Separation of cavity variable from others 

It is sufficient to replace R^'^ by the same quantity with the m-th term omitted: define 

Rh'^ := i y ala^ = R''' - ia^a^ = R''' + O f i V 
and similarly let 



Then we have the following relation, whose elementary proof is omitted: 
Lemma 2.8. For any v G [0, 1] and 6^ defined at it holds that: 

So = iQVrnfl ((^-L - r) (a^o-^ - r)) 

+ E 1^ E -'.0 m^/^o^ e--/- + O (iV-) . 

Thanks to a Lemma by Lattala and Guerra, we can now choose r in order to eliminate 
one of the terms in our overlap calculation Sq (at i; = with separated spins), as an 
immediate consequence of Lemma 12.61 

Lemma 2.9. For any choice of the parameters (3,'~fo, h > 0, the equation 

r = E [tanh^ {pZy^ + h)] (2.9) 
has a unique solution r G [0, 1], and we have 

i^m,o - r) {crl,al^ - r)) = 0. 

Let us take advantage of this relation, and try to write So in terms of r: going back to 
the expression in Lemma f2.8| we get 

*o = E I; E (^;:V.'-?)+o(^) 

+ 4 + (2.10) 
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with 



Al := S, V {Rf) , i E -^^0 - m e-"^ 



- is E ho {Ri') - " e 



%-Ki-k/N 



i&Cn 



where we have set 



Sk = e 



mi-k/N 



(2.11) 



We can bound now 5o in the following way: 



Lemma 2.10. Recall that T has been defined at Then there exists a constant k 

that depends on /? and q but is bounded for [3 bounded such that 



(2.12) 



Proof. Go back to relation ()2.1()|1 . and let us bound the terms A'l,A2, for k G Z'^. First 
of all, the estimation of A\ is controlled by 5*^. However, 



N 



N d 



ji=-N j^=-Nl=l 



n[(-l)''l(fc,^0)+iVl(,,=0) 



by an elementary argument on sums of geometric sequences. Hence, 

\^ I < ^ 

and thus, 



where Z(k) = number of components of k that are non-zero, (2.13) 



< 



kT 



(2.14) 



Furthermore, it is easily checked that 



\Al\<^ and |A||<./3V/^(|i?^'^f) J2^^^'^"{\^k 



(2.15) 



by applying inequality ()2.7|) in Lemma 12.51 Our claim is then proved easily by putting 
together pUjl and pT5|l . 

□ 

Remark 2.11. In the remainder of the article, k will stand for a positive constant de- 
pending on (3 and q and that is uniformly bounded in the range of our parameter (3; we 
will allow K, to change from line to line 
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Notice that, for our result on the fluctuations of Z^, we will need an improved bound 
on Al. This can be achieved under the following additional condition: 

Hypothesis 2.12. Going back to the decomposition \2.i3^) . we assume that there exists 
an integer d > d/2 such that, for every k E such that •jk 7^ 0, the number Z{k) of 
components of k that are non-zero satisfies Z{k) > d. 

Then Lemma f2.1UI can be enhanced in the following way: 

Corollary 2.13. Assume q satisfies Hupothesis \2.1\ and \2.1^ Then 

Uezd J/ 

Proof. This is trivially checked by going through the computations of Lemma f2. lUI again, 
and taking into account ()2.1Hj) . 

□ 

2.2.2 Task (ii). Difference between the overlaps at = and v = 1 

We are ready to state and prove the result which completes task (ii). 

Lemma 2.14. There exists a constant k, that depends on (3 and q but is bounded for (3 
hounded, such that 

\5i - (5o| < +Kl3^V J2 Ik^ (l^fcT) • (2-16) 
Proof. We can first write, using Lemma f2. 51 



i&Cn 



>1,2|2 



.kei 



(2.17) 



Similarly we can obtain 

70 \u {{R^'^ - r) {alal - r)) - u,, {{R^'^ - r) {a^ai - r)) | 



'^^/'f\R'''\' 



(2.18) 



We now get, by putting together ()2.17|) and (j2.18p and applying Jensen's inequality, that 



\6i-6o\=K(3^ 
which proves the lemma. 



□ 
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2.2.3 Self-averaging overlap limit 



We are now in a position to estimate 6i. We show tliat for small this expected total 
overlap, which is recentered using the value r, converges to at the speed as long as 
g is a continuous function. 

Proposition 2.15. Let (3 > and let r = r {(3) he the solution of \2.y\) . Let n he the 

constant defined in Lemma \2.8i and Lemma \2. 14\ i-e. n is a constant that depends on (3 
and q hut is hounded for [3 hounded. Assume that q satisfies Hupothesis \2.1l and that [3 is 
so small that 2k(3'^T < 1. In that case, we have, for N large enough, with R^f defined hy 
relation \2.'^) 

rV 



2k(3'^V)N 



Proof. We have, using (jTT^ and (jTTK|l 

< (5i = z/ 



O (i) . (2,19) 

where we recall that 0{N~^) is a function that tends to zero as fast as A^~^ and that 
this convergence holds uniformly in all parameters. Moreover, since k is bounded for (3 
bounded, for (3 sufficiently small we can make 2K/3^r smaller than 1. The result of the 
proposition follows. 

□ 

Corollary 2.16. Under the same assumptions as in Proposition \2. 1^ hut assuming addi- 
tionally that condition \2.12i holds true, then the conclusion of Proposition Wl^ holds with 
N replaced hy N'^. 

Proof. This follows trivially from the proof of Proposition 12. 151 if we modify the argument 
in order to take into account Corollary 12.131 

□ 



2.3 Consequence for the partition function 

We can now apply the previous computations in order to get the simple limit of Pn{i3), 
which recovers, in the high temperature region, the results contained in 

Theorem 2.17. Under the hypotheses o f Proposition \2A^ we have 
where the constant C depends on h, q, and (3, and is hounded for j3 G [0, (3q\ . 



< 



N ' 



(2.20) 
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Proof. Recall from Proposition 12.21 that 



1,2 I 2\ (i 



-,,.[\B^-' + r\')+^^^ (2,21 



The first term on the right-hand side can be handled easily: indeed, we have 

..2 , 



- y y 

J2d ^ ^ 

(«,i)GCjv feeZ" 

)-yy 

cy J\T2d ^—^ ^—^ 

^ 7o+ yik\s^ 



^^^iTTk-i/N ^-mk-j/N 



where has been defined at ()2.1ip . Thus, 



y oUi- 



J) 



(*j)eC]v 

and Proposition 12 . 1 51 then easily yields, for P < Pq, 



< 



K 



dpPN 



<^ + 7or/5|.(i?f)| 



2 2 

Furthermore, it can be shown, along the same lines as in jU], that for j3 < /5o, we have 



(2.22) 



\u{Rl'')\<K{P,q)/N, 



where K (/?, q) is bounded for P bounded. This inequality and (j2.22p now imply 



9f3PN - — (1 - r ) 



< — 

- N 



where k depends only on P,h,q and is bounded for P G [0,/9o]- The theorem follows 
by integrating Oppj^i, and using trivial calculations and known facts about the function 
SK. □ 



The final result we present in this section shows that while the complete structure of q 
does not seem to effect the limiting behavior of the partition function beyond the average 
value 7o of g^, the speed of convergence towards this value may depend heavily on the 
behavior of q. We show that the speed can be increased to the order A^~'^ as long as the 
hypotheses of Corollarv 12.161 hold. 

Corollary 2.18. Under the hypotheses of Corollarv \2.1b\ we have 



PN {P)-SKUj'p,h 



< 



where the constant C depends on h, q, and P, and is hounded for P G [0,/5o] 
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Proof. In the proof of Theorem 12.171 some estimates are aheady of order A^""^. For 
the others, we may use Corollary 12.161 instead of Proposition 12.151 in all its occurrences. 
This improves all estimates that were originally of order A^~^ to the order A^""^, with 
the exception of the estimation of the first term on the right-hand side of (j2.2ip . But 
here again, one can easily check that this term is of order A^~'^ under the conditions of 
Corollary EH 

□ 



3 Fluctuations of the free energy 

In this section, we will turn to our main aim, that is the central limit theorem governing 
the fluctuations of Zj^f. More specifically, we will get the following: 

Theorem 3.1. For (3 small enough, t G [0, 1], if q satisfies Hvpothesis \2.1\ and \2.1}A then 

(£) - hm [p^^^) _ p^^^ ^ 

where Y is a centered Gaussian random variable with variance t, and r is given by 

f : — , where f = Var [log(cosh(/3^7or2; + h))] (3.1) 



2 

with r defined by \2.9\} . and a standard Gaussian random variable z. 

This kind of result is usually obtained by letting all the interactions between spins 
tend to at once, and this procedure can be somewhat simplified by considering the 
computations from a stochastic calculus point of view. This leads us to consider a new 
path t I— i> Hj^^ti^cr) defined for t G [0, 1] by 

-H^A^) = ^ B,,,{t)qN{^ - j)a,a, + Pr'/^ ^ X,{t)cr, + ha„ (3.2) 

where f = 7or and {Bi j]{i,j) G C^} is again a collection of independent standard 
Brownian motions, and {Xf, z G Cat} is a family of independent reversed time Brownian 
motion, which can be seen as the solution to some stochastic differential equations of the 
form ()2.4|) . for a family {rii;i G Cn} (resp. {Wi;i G Cn}) of standard Gaussian random 
variables (resp. of independent Brownian motions). Here again, we will assume that all 
these objects are defined on the same probability space {Q,!F,P). Notice that, in order 
to spare notations, we have called this modified Hamiltonian Hp^ t again, like in Section |21 
hoping that this won't lead to any confusion. We will also denote by Zp^(t) and Pt{f) the 
partition function and the Gibbs average associated with HN^f Eventually, for t G [0, 1], 
we define ^ 

Pf3,h,t = (1 - rf + log(2) + E [log(cosh(/5y^^ + h))] , 

and we observe that Ppxi = Once this path has been defined, the strategy which 

leads to Theorem 13. II can be summarized as follows: apply Ito's formula in order to: 
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1. Compute the variations of t i— > e ^N,ti(r)_ 

2. Get an equation for the evolution of ¥^{1) = iV'^/^[p7v(/9) — Pf3,h,t\ 

3. Find a hmit for E[e*"'*^^(*)] for any ueR. 
We will detail now this global strategy. 

3.1 Preliminary computations 

Let us first study the dynamics of 1 1— e"^'^''*-'^^: 
Proposition 3.2. For t G [0, 1] and cr G S^v, we have 



exp ( a,{Pr'/\ + h)) + F^t) + F'^{t), 



with 



and 



idCN ° ^ 



Proof. The exponential function being a function with a nicely controlled growth, we 
can apply Ito's formula and we obtain 

Jo ' 2 Jo 

where (M)^ denotes the quadratic variation process of a semi- martingale M. Now, we 
can evaluate the quantity {H]\f^.)s, since all of the Bij, Wi are independent and using the 
fact that any finite variation process have a null quadratic variation. We get 

= 1^ E qU^ - J)t + P'fN% 

(ij)6Civ 

from which our claim is easily shown. 

□ 
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We will turn now to the fluctuations of Pn{P)- Namely set, for t G [0, 1], 



and define also the function $ : M ^ R by 

= log(cosh(/3f^/^a; + h)). 
Then the semi-martingale Yn can be decomposed in the following way: 
Proposition 3.3. Recall that T = 'Ylik^i'^'^k- Then, fort G [0, 1], Yj^{t) satisfies: 

l<2 

where the random variable and the processes Mi^]^ and Vk^N are defined by: 



(3.3) 



M2,N{t) 

^l,7v(t) 
V2,N{t) 



jj ft 

9N{i-j) / Ps{criaj)dBij{s) 



Nd/2 
Nd/2 



1 — s 



-ds 



J\fd/2 J\j3d/2 
jSfd/2p2 



E - J) - loN'/' 

■t 



0H 
4 



4 



fcezd J^ 



Proof. Notice that Ziq{t) is almost surely a strictly positive random variable. Thus, Ito's 
formula can be applied to log(Zjv(t)), and we obtain 



log(Zjv(t)) = log(Z^(0)) + 



*rfZjv(s) 1 d{ZN) 
Zn{s) 2 



^N^^) 



(3.4) 



The first two terms in the right hand side of ()3.4j) are easily computed. Indeed, it is easily 
checked that 



log(Z;v(0)) = NHog2+Y^ log[cosh(/?fi/2^i + h)]. 



(3.5) 
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Furthermore, invoking the fact that Zj\f{t) = XIo-gEjv ^ ^N.tio-] g^j^^j Proposition 13.21 we 
have 

Jo Zr,[s) N^l\^^c^ Jo ,^Jo 

Thus, putting together ()3.4p . (j3.5|) and ()3.6p . we have obtained that 

= [/^ + Mi,^(t) + M2,Nit) - V^^N{t) + V2,Nit) - / (3.7) 

2N"-i^ Jo ^Tvl'^J 



where 



(«j)eC]v 



Let us compute now the term (i(Zjv)s/Z^(s): according to Proposition 13.21 Z^ii) is a 
continuous semi-martingale, whose martingale part is 



MN{t) 



Hence, 



^Af(s) Jo ^Ifi^) 

i&CM-'^ {i,j)eCN 
Recall now that can be decomposed into g^(x) = X^fcez^ Thus, 

^^n.. Jo 



In order to simplify this last expression, we will use the elementary identity 

2 

E 
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valid for any family of complex numbers, and applied here to Recalling 
furthermore the definition ()2.3|) of -R^'^, we end up with 



^AtI-^) 



— 2^7fc / Ps [\Rk I jds 



Let us introduce now artificially the quantity -Rq^, by writing (i?^'^)^ = (_Rq'^)^- 
r^. This gives 



^Af('^) 



+ 2/3'iV'^f / p,{R'''')ds. (3.9) 



l.2^ 



Similarly to [Hj, let us notice that, since Xi{s) ~ J\f{0,l — s), a simple Gaussian 
integration by parts yields 



E 



Ps[cri 



1 - s 



E [9x,(.)P.(a.)] = Pr'/'E [p,(l - a.V,^)] . (3.10) 



This elementary consideration will induce us to add and subtract SiPs(l ~ '^I'^u to 
the expression (j3.7|) . By plugging moreover (j3.9|) . we get 

Y^{t) = Um + M.^^it) + M2,^(t) - {V^,N{t) - V2,N{t)) + V,,Nit), (3.11) 



where 



Nd/2 



ft q2 

y / ps{l- (yla'^i)ds+ : ^ V qy,i-j)t 



4:Nd/2 



+ 



f / p,{R'^')ds- 



l-r) 



This last expression can be simplified a little to give 

(«j)eCjv 



4 



that is V3,Ar(t) = V3^7v(^)- By reporting this equality in (|3.1ip . the proof is now complete. 

□ 
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The last preliminary result we will need in order to establish our CLT is a self-averaging 
result for i?^;^ under the measure pt at a fixed value of the parameter t G [0, 1]. 

Proposition 3.4. Under the hypotheses of Corollarv \2.16{ let t G [0, 1]. Then 

1 



(3.12) 



for d > d/2, uniformly in t. 

Proof. This is an easy elaboration of the computations leading to Proposition 12.151 by 
considerering the path t> G [0,t] HN^t,v{o') defined by 



with obvious notations. 



□ 



3.2 Proof of Theorem [O 

For an arbitrary m G M, we will try to control $Ar^^j(t) = E[e*"-*^^'^*)]. To this purpose, we 
will apply Ito's formula to the complex valued Cl function x i— e*"^. We obtain, for any 
tG [0,1], 



where 



and 



^iuYN{t) 



m=2 



iul3 



D5,N{t) 



iu[3f 



1/2 



Ps{cri)dWi 



lU 



1 — s 



fjd/2 



(ij)6Cjv 



9M^-JJ -70 



fjd/2 



pjd/2p2 
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The terms D-i^N{t) and Ds^N{t)a.ie the ones associated with the quadratic variation process 
of Y]y(t), that is: 







n U P r ^ /■ j„yjy(s)^ /'^1^2\ U P f f iuYn(s)^ /'d1,2\j, 

Ds,N[t) = — —^2^ / e ^^V4aiajds = — / e ^^>ps{R')ds. 

i&Cn ° ° 

Let us find some estimates for the expected value of all the terms we have obtained in 
this decomposition. First of all, the usual central limit theorem for IID random variables 
yields 



(3.13) 



where f is defined at ()3.H) . Furthermore, the terms -D2,Af(^) and D^^N{t) are obviously of 
zero mean. 

In order to control D^^^it), let us perform again the Gausssian integration by parts 
()3.1()|1 . which can be read here as 



]Sfd/2 
2/9^:1/2 



u j3r 



J\fd/2 



/ E[9x,(.)(F^(s))e»^-Wp,(a,)]ds-E[D5,^(t)]. 



Furthermore, 



and therefore, we obtain 



Jo 

= u^(3''^y ^^^^[s)ds + o(^-^^ 

owing to Proposition 13.41 



(3.14) 



Let us turn now to the estimation of DQ^N{t): notice that, under Hvpothesis 12.11 and 
12.121 it is easily checked that 



fT2d XI 

(i,i)6Cjv 



.^ 70 

^-j)-y 
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which together with a direct apphcation of Proposition 13.41 shows that 

where e = d — d/2. 

As far as Dj^N{t) is concerned, notice that by relation ()3.8|) . we have 



4 

and thanks to proposition 13.41 we get 



/ e^"^-Wp.(|i?''T)^s + ^^ / e»^-Wrfs, (3.15) 
Jo AN"^ Jo 



DMt) = -"^^ fe--^^^^ds + 0(-^^\ , (3.16) 



and thus 



Eventually, in a similar way we have 

E [DsMt)] = -"^^^ f <^N,u{s)ds + O f-^ 
^ Jo \N"- 

Putting together the previous estimates on E[Di^N(t)], . . . , E[Ds^Nit)], we have finally: 

4 Jo 

with 

K 



RN,u(t) 

which ends the proof by a Gronwall type argument. 



< 



□ 
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